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Abstract 


Using  six-vector  formalism,  the  polarisability  problem  of  a layered  bi-anisotropic  ellipsoid  is 
solved.  The  polarisability  six-dyadic  is  explicitly  calculated  for  the  case  of  a two-component 
ellipsoid  where  the  core  is  fully  bi-anisotropic  but  the  shell  and  the  environment  are 
bi-isotropic  at  most.  It  is  reasonably  straightforward  to  include  more  bi-isotropic  layers  on 
the  ellipsoid.  The  limitation  of  the  analysis  is  that  all  ellipsoidal  boundaries  for  the 
composite  structure  have  to  be  confocal. 


1.  Introduction 

In  this  presentation  we  treat  and  analyse  the  calculation  of  the  polarisability  of  a small  scatterer. 
The  polarisability  is  defined  as  the  linear  relation  between  the  dipole  moment  and  the  uniform 
incident  field  that  induces  this  moment.  For  anisotropic  or  non-spherical  scatterers,  the  polar- 
isability is  generally  a dyadic:  p = E ■ E,  where  p is  the  dipole  moment  and  E is  the  incident 
electric  field.  In  the  present  paper  we  will  concentrate  on  the  case  of  bi-anisotropic  scatterers 
in  which  case  the  polarisability  is  a matrix  relation  between  the  induced  electric  and  magnetic 
dipole  moments  and  the  incident  electric  and  magnetic  fields.  The  concept  of  polarisability  is 
important  in  the  low-frequency  applications  [1]:  the  size  of  the  scatterer  has  to  be  small  in 
comparison  of  the  wavelength  of  the  operating  field. 

In  the  bi-anisotropic  case  we  will  use  the  six- vector  formalism  [2].  The  electric  and  magnetic 
vector  quantities  are  collected  into  a six-vector,  and  the  relation  between  two  six-vectors  is  a 
six-dyadic.  Thus,  for  example,  the  polarisability  is  a six-dyadic  (consisting  of  four  ordinary 
three-dyadics,  or  36  scalar  parameters):  p = A -e,  where  the  field  and  dipole  moment  six-vectors 
are 

p=(pI)  and  e=(n)  « 

In  this  presentation,  the  particular  emphasis  is  on  inhomogeneous  scatterers;  especially  lay- 
ered ellipsoidal  scatterers.  We  will  show  how  the  polarisability  six- matrix  of  a layered  ellipsoidal 
scatterer  can  be  calculated  with  the  six-matrix  formalism:  the  equations  are  shown  for  the  case 
of  core-plus-shell  ellipsoid  but  the  method  works  for  any  number  of  ellipsoidal  layers.  There 
is,  however,  one  limitation:  all  the  ellipsoidal  boundaries  in  the  structure  have  to  be  confocal. 
In  other  words,  the  foci  of  the  various  ellipsoids  have  to  coincide.  The  results  allow  the  layers 
(and  the  environment)  to  be  bi-isotropic;  in  addition,  the  core  of  the  layered  ellipsoid  can  be 
arbitrarily  bi-anisotropic. 
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Figure  1:  Layered  ellipsoid  with  the  core  and  shell 


The  results  form  a considerable  generalisation  of  the  previous  studies  where  the  corresponding 
polarisability  problem  for  layered  ellipsoidal  dielectric  scatterers  [3]  and  for  the  two-layer  chiral 
sphere  [4]  was  solved. 


2.  Polarisability  Six-Dyadic  for  the  Core-Plus-Shell  Ellipsoid 

Consider  the  scatterer  illustrated  in  Figure  1 where  an  ellipsoid  is  located  in  bi-isotropic  envi- 
ronment. The  ellipsoid  consists  of  a bi-anisotropic  core  and  a bi-isotropic  layer  (shell).  It  has  to 
be  remembered  that  the  ellipsoidal  boundaries  have  to  be  confocal,  which  roughly  means  that 
the  core  is  less  ’’rounded”  than  the  whole  scatterer.  The  material  parameters  in  the  problem 
are  contained  in  the  three  material  six-dyadics 
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How  to  calculate  the  six-dyadic  polarisability  for  this  ellipsoid?  Using  the  three-vector  anal- 
ysis as  was  done  in  [4]  for  the  layered  chiral  sphere  would  lead  to  a large  set  of  coupled  vector 
equations  from  which  the  polarisability  would  be  extremely  tedious  to  solve.  Here  the  six-vector 
analysis  shows  its  power,  because  the  analysis  remains  formally  on  the  same  level  of  complexity 
as  in  the  three- vector  case,  and  we  can  follow  the  steps  in  the  layered  dielectric  studies. 

The  result  is  that  the  polarisability  six-dyadic 


A = 


aem 

amm 


(3) 


of  a layered  bi-anisotropic  ellipsoid  can  be  expressed  in  a compact  form: 

A = VMe  • {K  - [(Mi-Me)  • h + Me]  + «7(M2-M1)  • U}-1  • [K  • (Mi  — Me)  + w(M2-Mi)]  (4) 


c 


with  the  following  definitions: 


K = (M2-M0  • (l2-w\-i)  ■ MJ-1  + I 


and  the  unit  six-dyadic  is  defined  obviously  by 
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The  depolarisation  factors  of  the  two  ellipsoidal  boundaries  axe  contained  in  the  six-dyadics 

£)  and  L2  = (*»  £) 

V»  Li)  \0  L2J 

and  w is  fraction  of  the  core  from  the  total  volume  V of  the  particle. 

The  present  analysis  can  be  generalized  to  ellipsoids  having  more  layers  than  the  core  and 
the  shell,  as  in  [3,  5]. 


3.  Special  Cases 

3.1  Homogeneous  sphere 

It  is  important  to  check  the  result  with  special  cases  known  earlierly.  A trivial  test  could  be  to 
have  a homogeneous  isotropic  sphere  in  isotropic  environment,  for  which  the  scalar  polarisability 
is 


q = 3Vee 


d -ce 

d -f  2ee 


It  is  quite  easy  to  see  that  (4)  passes  this  test. 


(5) 


3.2  Layered  dielectric  ellipsoid 


The  second  special  case  is  the  layered  dielectric  (or  magnetic)  ellipsoid  (see  [3]):  Ce  = Ci  — Ce  = 
Ci  — 0)  C2  ==  C2  = 0,  C2  = 62 /,  and  /Z2  — M2 /•  In  this  case  the  polarisability  six-dyadic  reduces  to 

A /«ee  =0  \ 

\ 0 amm ) 

where  the  electric  polarisabihty  is 

= _ ye  Y'  (Cl  - d=)(d  + N2i(e2  - ei))  + n?(ei  + Aii(ee  - d))(e2  - d.) 

“ 6 ^ (ee  + Nu(e  1 - 6e))(ei  + AT2i(e2  - d))  + wNu(l  - JVu)(d  - ee)(e2  - d)  "iUi 

(6) 

and  the  magnetic  polarisability  SmTn  is,  mutatis  mutandis,  the  same  as  See.  Nu  and  N2{  are 
the  depolarisation  factors  of  the  ellipsoids: 

-^1  = NlXUXUX  + NlyUylly  + N 3nd  Is2  = N2xUXUX  + N2yUy\ly  + N2z  uz  uz. 


3.3  Other  special  cases 

Among  the  other  possible  special  cases  to  test  the  result,  the  following  ones  can  be  considered: 
homogeneous  chiral  sphere  and  layered  chiral  sphere.  The  result  (4)  passes  also  these  tests,  but 
the  resulting  lengthy  expressions  are  not  shown  here. 

4.  Maxwell  Garnett  Mixing  Formula 

With  (4)  we  can  build  the  Maxwell  Garnett  (Clausius-Mossotti)  mixing  formula  for  aligned 
layered  ellipsoids: 

= Me  + nMe  • ^Me  — nA  • Li^  -A 
where  n is  the  number  density  of  the  layered  ellipsoidal  inclusions. 


(7) 
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Figure  2:  Effective  relative  permittivity  (left)  and  chirality  (right)  of  the  mixture. 


Numerical  Example 

Let  us  assume  a mixture,  which  consists  of  layered  chiral  spheres  in  an  isotropic  background. 
The  cores  in  the  spheres  are  bi-isotropic:  ?2  = 2.0eo I,  P2  — 1.5po I>  and  «2  = —1-0,  so  that 
£2  = jy/UQtoI  and  C2  = The  shells  are  also  bi-isotropic:  = 2.0eoI,  /+  = 1.5hqI, 

but  /«i  = +1.0.  In  other  words:  the  materials  in  the  layered  spheres  are  almost  the  same,  but 
the  cores  are  right-handed  and  the  shells  are  left-handed.  For  the  background  medium:  ee  — to 
and  ne  = fifl. 

It  is  a quite  simple  task  to  write  a short  Matlab-code  and  use  (4)  and  (7),  to  calculate  the 
effective  material  properties  of  the  given  mixture.  In  Figure  2 is  shown  the  effective  relative 
permittivity  and  chirality  of  the  mixture.  The  volume  fraction  / is  defined  as  f — nV , where  n 
is  the  number  density  of  the  inclusions  in  the  mixture  and  V is  the  total  volume  of  an  inclusion. 

Figure  2 shows  that  although  both  components  of  the  spheres  have  the  same  electric  per- 
mittivity, the  mixture  permittivity  eeff  is  not  the  same  in  the  case  / = 1 (no  background, 
everything  just  inclusions).  The  strange  effect  on  permittivity  is  caused  by  the  magnetoelectric 
coupling. 
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